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ABSTRACT 

Period-colour (PC) and amplitude-colour (AC) relations at maximum, mean and 
minimum light are constructed from a large grid of full amplitude hydrodynamic mod- 
els of Cepheids with a composition appropriate for the SMC (Small Magellanic Cloud). 
We compare these theoretical relations with those from observations. The theoretical 
relations are in general good agreement with their observational counterparts though 
there exist some discrepancy for short period (log[P] < 1) Cepheids. We outline a 
physical mechanism which can, in principle, be one factor to explain the observed 
PC/ AC relations for the long and short period Cepheids in the Galaxy, LMC and 
SMC. Our explanation relies on the hydrogen ionization front-photosphere interaction 
and the way this interaction changes with pulsation period, pulsation phase and metal- 
licity. Since the PC relation is connected with the period-luminosity (PL) relation, it 
is postulated that such a mechanism can also explain the observed properties of the 
PL relation in these three galaxies. 
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1 INTRODUCTION 

The papers in this series are devoted to a study of the period- 
colour (PC) and amplitude-colour (AC) relations for classi- 
cal Cepheid variables, the physics behind the connections 
of PC and AC relations, and the implication for Cepheid 
period-luminosity (PL) relations. The connection betwee n 
the PC and AC relations was derived bv lSimon et alj l|l99&t ). 
who applied the Stefan-Boltzmann law at the optical maxi- 
mum and minimum light together with the f act that radial 
variations are small in the optical (|Coxlll980l ). Specifically: 



log Tmax log Tjriin (Vrn 



V m 



(1) 



where T max and T min are the effective temperature at the 
maximum and minimum light, respectively. If T max , and 
hence the colour at maximum light, is independent of pe- 
riod, then equation (1) implies there is a relation between the 
temperature at minimum lig ht and the optical amplitudes, 
and vice versa. Code ( 1947) found that at the maximum 
light, the spectral type of Galactic Cepheids is indepen- 
dent of pulsating period. Convincing observational support 
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for a flat (i.e. zero s lope) maximum lig h t PC relation has 
been d ocumented by ISimon et alj (1 19931) . Kanbur fc Ngeowl 



|2004 here after Paper I) jKanbur et all l|2004 hereafter Pa 
per II) and lKanbur fc Ngeowl (|2006l . hereafter Paper III) for 
the Galactic (with log[P] > 0.8, where P is pulsation pe- 
riod in days) and LMC (with log[P] > 1.0) Cepheids. These 
authors imply that the flatness of the Cepheid PC relation 
at maximum is due to the interaction between the hydro- 
gen ionization front (HIF) and the photospher e (defined as 
optic al depth r = 2/3; see Paper II, III and ISimon et all 
Il993l) . However, Paper I found that the SMC PC relation, 
in (V — I) colour, at maximum light is not flat, but in- 
stead suffers a "flattening" of the slope from 0.396 ± 0.039 
to 0.207 ± 0.071 for short and long period (log[P] > 1.0) 
Cepheids, respectively. This change of slope is marginally 
significant using the P-test described in Paper I. One moti- 
vation for this paper is to investigate the physics behind PC 
relations at maximum light and in particular why Galactic 
and LMC Cepheids display flat maximum light PC relations 
and SMC Cepheids do not. 

Since the PC and PL relations are really projections of 
the period-luminosity-colour (PLC) relation, PC and PL re- 
lations are closely connected: changes in one are generally 
reflected in changes in the other. Hence one way to look for 
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changes in Cepheid PL relations is to look for changes in PC 
and AC relations. Another important driving force behind 
this series of papers is to investigate the possible physics be- 
hind recent results suggesting that the LMC (Large Magel- 
lanic Cloud) PL (and PC) relation is non-linear: in the sense 
that current LMC Cepheid data are statistically consistent 
with two PL (and PC) relations of significantly differing 
slope with a break at/around a period o f 10 days. The obser- 



vation al evidence for th is is presented i nlTammann fc Reindl 
2002f) . Pa per I fe III. INgeow et al.l (|2005h . ISandage et"al 



20041 ) and INgeow fc Kanburl (|2006al . hereafter Paper IV) 



It is also interesting to note, as has been pointed out in 
previous work, that the "break" period of log(P) ~ 1.0 is 
also the location of the well known Cepheid pulsation reso- 
nance. Paper IV use the OGLE (Optic al Gravitational Lens- 
ing Experiment, lUdalski et aT]|l999bl ) LMC data and phase 
all Cepheids to a common epoch and then plot both PC and 
PL relations as a function of phase. This approach presents 
clear and compelling evidence of a break at a phase around 
0.8 and also attest to the dynamic nature of the PL rela- 
tion. The PC relation follows changes in the PL relation 
very closely. However the standard PL relation used in the 
literature is at mean light: the average over phase of the 
multiphase relations presented in Paper IV. Because the PL 
relation is indeed linear at certain phases (e.g., at maximum 
light), the net result is that the strength of the non-linearity 
is diluted at mean light. We emphasize, though, that a bat- 
tery of parametric and non-parametric statistical tests car- 
ried out on the LMC PL relation at mean light in opti- 
cal bands confirm the existence of the non-linearity at high 
confidence levels - greater than 99%. These tests include 
the F-test, robust methods like Tukey's bi- weight func- 
tion and non-parametric techniques (Paper I; INgeow et al.l 
120051 ). In contrast the Galactic and Small Magellanic Cloud 
(SMC) PL an d PC relations are f o und to be linear with 
current data (|Udalski et al.l [l9~99al: iTammann et al.l 120031 ; 
iKanbur fc Ngeowll2004l ; INgeow fc Kanburl 120041 ). Since the 
Galaxy, LMC and SMC have different metallicity, it has been 
conjectured that metallicity is the key to understanding this 
behavior. 

We now briefly discuss some criticisms of results imply- 
ing a non-linear LMC PL relation. One of these is that the 
OGLE LMC sample lacks sufficient long period Cepheids. 
Despite the fact that the F-test takes int o account both the 
number and nature of the observations, INgeow fc Kanbur] 
(2006b) augmented the OGLE sample with data taken 
from, for example. ICaldwell fc Lanevl |l99lf) and lSebo et alj 
120112 ). They conclude on the basis of their rigorous statisti- 
cal test (the F-test) that even with this augmented LMC 
Cepheid data, the LMC PL relation is still found to be 
non-linear in the sense described above (even though the 
difference in the slopes of the long and short period PL 
relations may be small). Another criticism is that extinc- 
tion errors are influencing the result even though only pub- 
lished extinction values are used. Well established observa- 
tions and theory have shown that the Galactic and LMC 
Cepheid PC relations are flat at maximum light, at least for 
the lon ger period (log[PJ > 1) Cepheids (see Paper I, II, 
III and lSimon et al.||l993T ). If the extinction values required 
to make the mean light LMC PC and PL relations linear 
are adopted, then the LMC PC relation at maximum light 
is no longer flat. In this scenario, LMC Cepheids get hot- 



ter at maximum light as the period increases, which is in 
distinct disagre ement with the th eory and observations at 
maximum light l|Simon et al.l ll993) . The referee has pointed 
out that if the amplitude increases sufficiently fast with pe- 
riod, then it could produce hotter temperatures at maximum 
light with increasing period. But note that the amplitude in- 
creases fairly rapidly after 10 days and this is precisely when 
the observed PC relation is flattest. Further if the LMC 
Cepheids are such that they got hotter at maximum light as 
the period increases this is different to Galactic Cepheids at 
maximum light. We therefore believe that extinction errors 
are not the cause of the observed non-linear LMC PL rela- 
tion. More detailed discuss i ons r e garding extinction e rrors 



are given in ISandage et ail l|2004) . INgeow et all i|2005l ) and 



INgeow fc Kanburl feooebn T and will not be repeated here 



Since the Cepheid PL relation is a cornerstone of mod- 
ern astrophysics, a proper understanding of its properties 
and the physics behind it is crucial. For example, it is of fun- 
damental importance in establishing an extra-galactic dis- 
tance scale that is independent of CMB studies. An appro- 
priate refinement of the Cepheid PL distance scale, can pro- 
vide a value of Hubble's constant accurate to less than 3-4%. 
Though the current Cepheid base d estimate of Ho is accu- 
rate to 10% (|Freedman et al.l l2001). a more accurate (< 5%) 
direct measurement of Ho via a Cepheid distance scale is 
still very relevant because it will help to break the degen- 
eracy between fimattor and Hp present from WMAP CMB 
estimates (see , for example, Tegmark et al.ll2004l ; iFreedmar] 



l2005l : lHul2005l : ISper g el eta! 



20071 ) . CMB measurements can 



only estimate f2 ma tterflo an d estimate Ho only if a sim- 
ple fla t ACDM model is assumed. Table 2 of ISpergel et al.l 
l|2007l ) points to the fact that an independent estimate of 
Ho accurate to 1% will result in a reduction of the 65% con- 
fidence interval on f2 ma ttor by almost a factor of two over 
that with WMAP alone. This is a greater red uction than if 
other data such as 2dFGRS are used (table 5 of lSpergel et al.l 
l2007h . With the possibility of Gaia and SIM satellites pro- 
ducing more accurate parallaxes for many Cepheids a con- 
siderably more accurate Cepheid zero point is well within 
reach. Better refinement of the slope of PL relation will be- 
come vital. Equally as important, modeling the Cepheid PL 
relation can yield important information about the mass- 
luminosity (ML) relation obeyed by Cepheids. Since these 
ML relations are sensitive functions of input physics, such 
as the amount of convective overshoot, such modeling is also 
important for theories of stellar evolution and pulsation. 

A working theoretical hypothesis to explain the non- 
linear LMC PL and PC relation is the interaction of the 
HIF with the stellar photosphere at certain phases. Paper II 
& III published theoretical pulsation models of Galactic and 
LMC Cepheids, respectively, to model this behavior. In these 
papers we found that the interaction of the HIF and stellar 
photosphere may play an important role in explaining the 
observed PC and PL behavior. We briefly summarize this 
theory in what follows (see Paper II & III for more details) . 
In a pulsating Cepheid, the HIF and stellar photosphere in- 
dependently move in and out in the mass distribution as 
the star pulsates. In certain circumstances, the stellar pho- 
tosphere can be located at the base of the HIF. In this case, 
when the the density is low the temperature of the stel- 
lar photosphere and hence the effective temperature of the 
Cepheid at that phase experiences a much reduced depen- 
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Table 1. Input parameters for SMC Cepheid models with periods 
obtained from a linear analysis. The periods, Po and Pi, are re- 
ferred to the fundamental and first overtone periods, respectively. 
Similarly for the growth rate, rj. Both of the mass and luminosity 
are in Solar units, the temperature is in K and the period is in 
days. 
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dence on pulsation period. Thus the PC relation is flatter 
in this situation. Because the PL relation is intimately con- 
nected to changes in the PC relation, the PL relation slope 
will also change. Because this HIF-photosphere engagement 
happens suddenly, the changes in the PC and PL relation 
are sharp. The phase and period at which such a low density 
HIF-photosphere interaction can occur depends on the ML 
relation since this controls the location of the HIF in th e 
mass distribution (|Kanburlll995l ; iKanbur fc Phillips! 1 199rj ). 
The connection to metallicity occurs because the ML rela- 
tion, taken from stellar evolution calculations, is indeed a 
function of metallicity. 

In this paper we extend the work of Paper II & III by 
constructing a large grid of full amplitude hydrodynamic 
models of Cepheids with an SMC metallicity and study how 
the interaction of the photosphere and HIF affects the PC 
(and hence the PL relation) in these models compared to 
previous work looking at Galactic (Paper II) and LMC (Pa- 
per HI) models. Section 2 outlines the methodology and pul- 
sation codes we used in this paper and Section 3 presents our 
results for the SMC models. The conclusion and discussion 
is given in Section 4. 



Table 2. Temperatures at maximum and minimum light from 
full-amplitude non-linear model calculations. The periods, lumi- 
nosity and temperature are in days, Lq and K, respectively. 



Lraax ^max L m { n T^vnin 



ML Relation from Bono et al. (2000) 



£9 aQK 


53454.42 


5988. 


56 


27021.41 


ouoo. 


31 

. ■)! 


44 130 


35895.77 


5826. 


.30 


19377.45 


5022 


77 


oo.yoi 


21806.85 


5480. 


18 


14111.82 


ouoi . 


9 A 


^i Q 


15662.90 


5380. 


53 


11737.37 


ouyu. 


. ol) 


99 ^fi.9 


11969.05 


5626. 


25 


8755.998 


OUOO. 


77 




7880.140 


5581. 


2.3 


6253.112 


oiijy. 


7 r » 


1 Q 91 Q 


5731.252 


5560. 


19 


4789.290 


OZDU. 


. >.);) 


1U.OUO 


3934.212 


5488. 


40 


3546.536 




-r. 

. DO 


Q £37 


3699.569 


5495. 


11 


3336.452 


Oo / o. 


7-1 


O.O 1 U-L 


3093.723 


5609. 


36 


2808.320 


5421 


22 




2555.649 


5607. 


.43 


2357.080 


044 ( 1 . 




7 1171 


2389.147 


5605. 


87 


2225.435 




on 


(\ 7999 


2234.535 


5625. 


98 


2089.173 


OOU i . 


R1 

Ol 


D.UOOO 


1943.932 


5642. 


11 


1839.759 


oooo. 


I {-I 
. 10 


5.4218 


1681.461 


5664. 


88 


1604.783 


5574. 


.01 




ML Relation from Chiosi (1989) 
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2 METHODS AND SMC MODELS 

The codes for computational and numerical me thods for do- 
ing th e pu lsation modelings ar e described in lYecko et al.l 
|l998t ) and lKollath et all (2002), and are exactly the same 
as used in Paper II & III. In brief, the input parameters 
to the pulsation codes include the mass (M), luminosity 
(L), effective temperature (T e ff) and chemical composition 
(X, Z). In this paper, the chemical composition is set to be 
(X,Z) = (0.70,0.004) to represent the SMC hydrogen and 
metallicity abundance (by mass). The input effective tem- 
peratures are chosen to ensure the models oscillate in the 
fundamental mode and located within the Cepheid instabil- 
ity strip. To be consistent with Paper II & III, we adopt the 
same ML relations as given below: 

(i) ML relation given in lBono etaLI (|2000h : 

log(L) = 0.90 + 3.35 log(M) + 1.36 log(F) - 0.34 \og(Z), 
= 3.35 log(M) + 0.996. (2) 

(ii) ML relation given in lChiosil l| 19891 ) for Z = 0.001: 

log(L) = 3.22 log(M) + 1.511. (3) 

The units for both M and L are in Solar units. To con- 
vert the temperatures from the models to the (V — I) 
colours, we use the BaSeL atmosphere database (jLeieunel 
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2002; Wcs tera et alj|2002j ) to construct a fit giving temper- 
ature and effective gravity (see Paper ff) as a function of 
(V — /) colour. The bolometric corrections (BC) are con- 
verted in a similar manner. Note that the effect of a micro- 
turbulence parameter that varies with phase has been dis- 
cussed in Paper III: its variation with phase does not in- 
fluence the results presented in this paper. To convert the 
observed colours to the temperatures appr opriate for the 
SMC data, we use the prescriptions given in lBeaulieu et al,l 
ll200ll) : 

log( S ) = 2.62- 1.21 log(P), 
log(T e //) = 3.91611 + 0.0055 log(g) - 0.2482(V - 7) , 
AT = log(T e// ) - 3.772, 
BC = -0.0324 + 2.01 AT - 0.0217 log(ff) 
-10.31(AT) 2 . 

More details of the methodology, the pulsation codes and 
the temperature-colour conversions can be found in Paper 
II & III, and will not be repeated in detail here. 



3 RESULTS 

The results for our SMC models are collectively summarized 
in Table Q]|3] Table [T]presents the input mass (and hence the 
luminosity) and effective temperature for our SMC models 
and the resulting periods from linear non-adiabatic calcula- 
tions. Table [T] displays the fact that the models constructed 
in this study occupied a wide range of periods with an unsta- 
ble and stable fundamental and first overtone mode respec- 
tively. Thus we avoid regions of the instability strip which 
are susceptible to first overtone or double mode pulsation. 
Table [2] gives the temperatures at the maximum and mini- 
mum luminosity for the models calculated in Table [T] Tabic 
[3] with identical layout as in table 3 of Paper II, presents 
more detailed information regarding photospheric temper- 
atures at mean luminosities during the ascending and de- 
scending branches of the light (or luminosity) curves (see 
Paper II & III for more details). 

The results from the SMC models were compared to 
the observations, also presented in Figure [TJ3] In these fig- 
ures, Fourier fits to the SMC Cepheid data, taken from Pa- 
per IV, were used to calculate the observed amplitudes and 
colours at maximum, mean, and minimum light using the 
definitions given in Paper I & III. Table [4] & [5] summarizes 
the empirical PC and AC relations at maximum, mean and 
minimum light for the SMC Cepheid data plotted in these 
figures, respectively. In Figure [TJ we compare the temper- 
atures given in Table & to the observed SMC Cepheid 
data. Similarly, Figure [2] displays the PC relations at various 
phases for the SMC models superimposed on the SMC data. 
From these figures, we found that the SMC PC relation at 
minimum light has the smallest scatter for all periods and 
does not suffer a statistically significant change of slope at 
log(P) ~ 1.0, i.e. it is linear over the entire period range. 
The PC relation at maximum light is significantly flatter for 
long (log[P] > 1.0) period than short period Cepheids and 
displays marginal evidence of a slope change at log(P) ~ 1.0. 
It also has the largest scatter over all phases. 

The SMC models do a reasonable job in matching 



the observations in the period-temperature and the period- 
colour plane, though the short period models at maximum 
light tend to be too cool compared to the observations (as 
in Figure [T] & [2]) . Figure [3] portrays the results of AC rela- 
tions from the models as compare d to the SMC data . Several 
models, especially those using the lBono et al.l (|2000h ML re- 
lation, display a smaller amplitude as compared to the data. 
This is seen in the Galactic (Paper II) and LMC (Paper III) 
models as well, however we believe this does not affect our 
results on the HIF-photosphere interactions (Paper III). To 
the best of our knowledge a comparison of models and the- 
ory on PC/ AC planes at different phases, as described here, 
has not been carried out in such detail before: further work 
will investigate in detail some of the discrepancies between 
models and observations described in the text. 

From the full amplitude SMC models, the temperature 
and the opacity profile can be plotted at a given phase of 
pulsation. Similar to Paper II & III, the locations of the HIF 
(sharp rise in the temperature profile) and photosphere can 
be identified from the temperature profile. Figure [4] displays 
the temperature and opacity profiles, with the locations of 
photosphere marked as filled circles, for a long period, a 10- 
days and a short period SMC model. From the temperature 
profile, the "distance" A between the HIF and the photo- 
sphere from the temperature profile can be calculated (see 
Paper II for the definition of A). A small A implies there 
is a HIF-photosphere interaction, and vice versa. Figure [S] 
presents the A as a function of pulsating period for the SMC 
models with the two ML relations used in this paper. In Pa- 
pers II & III, it is found that the distribution of A as a 
function of period is almost independent of the adopted ML 
relations. This is also seen in the SMC models as depicted in 
Figure^ Furthermore, we see, in contrast to the behavior for 
Galactic and LMC models, the distance is reasonably con- 
stant for all periods and at all phases. Figure 15 of Paper III 
shows that at minimum light this distance for the Galactic 
models increases with period, but for the LMC models this 
distance only increases for periods greater than 10 days. At 
maximum light, all models display a constant distance. The 
SMC models are different in that there is almost a constant 
distance between the HIF and stellar photosphere, indepen- 
dent of pulsation period, pulsation phase and ML relation. 

Figure [5] depicts the photospheric density as a func- 
tion of period for the SMC models. This figure needs to be 
compared with figure 16 of Paper III. The difference be- 
tween the SMC models and the Galactic/LMC models is 
in the photospheric density at minimum light. The Galac- 
tic models have a photospheric density at minimum light 
which is less than 10 _8 gcm -3 . Paper II also found that for 
these Galactic models, the photosphere and HIF are not en- 
gaged. For the LMC models at minimum light, Paper III 
found that the HIF and photosphere are nearly always en- 
gaged with some indication that there is a separation for 
periods greater than 10 days. We also see that the photo- 
spheric density for these longer period LMC models is less 
than 10 _8 gcm -3 . The SMC models always have an engaged 
photosphere and HIF at all phases. But crucially, the den- 
sity of the photosphere is higher than 10 _8 gcm -3 for mod- 
els with log(P) < 1.5. The higher densities mean that no 
sharp change in the temperature of the photosphere inde- 
pendent of the global stellar parameters can occur because 
our premise is that a sharp change can occur in the slope 
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Table 3. Temperatures at mean light from full-amplitude non-linear model calculations. See Paper II for the meanings of < L >, L 
T m ean and T™dan ■ The periods, luminosity and temperature are in days, Lq and K, respectively. 
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Figure 1. The plots o f log(T)-log(P) rel ation s for th e SMC data (small crosses) and models. The open and solid squares are for the 
models calculated with iBono et al. I feOOd) and IChiosil dl98Sh ML relat ions, respectively. The conversion of the (V — I) colours to the 
temperature are done using the equations given in lBeauIieu et al.l fl200lt) . Left panel: log(T)-log(P) relations at maximum and minimum 
light. Right panel: log(T)-log(P) relations at mean light for both of the ascending and descending means. 
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Figure 2. The SMC PC relations with the results from the models. The symbols are the same as in Figure[T] The temperatures of the 
models are converted to the (V — I) colour using the BaSeL database. Left panel: PC relations at maximum and minimum light. Right 
panel: PC relations at mean light for both of the ascending and descending means. 
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Figure 3. The SMC AC relations with the results from the models. The symbols are the same as in Figure [T] The bolometric light 
curves from models are converted to ^/-band light curves with the BC obtained from the BaSeL database. Left panel: AC relations at 
maximum and minimum light. Right panel: AC relations at mean light for both of the ascending and descending means. 



Table 4. The SMC period-colour relation in the form of (V — I) = alog(P) + b, and a is the dispersion of the relation. Long and short 
periods refer to Cepheids with log(P) > 1.0 and log(P) < 1.0, respectively. See Paper I & III for the definition of the phases. 
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Table 5. The SMC amplitude-colour relation in the form of (V — I) = aV am p + b, and a is the dispersion of the relation. Long and short 
periods refer to Cepheids with log(P) > 1.0 and log(P) < 1.0, respectively. See Paper I & III for the definition of the phases. 
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Figure 4. The temperature (top panels) and the opacity (k, bottom panels) profiles, plotted in terms of the internal mass distribution 
(log[l — M r /M], where M r is mass within radius r and M is the total mass), for a long period, a 10-days period and a short period SMC 
models. The dotted, solid and dashed curves are for the profiles at maximum, mean and minimum light, respectively. The filled circles 
denote the location of the photosphere at r = 2/3 for each phases. The mean light profiles at the ascending and descending branch are 
the solid curves that lie close to the profiles at maximum light (dotted curves) and minimum light (dashed curves), respectively. 



of the PC relation when the photosphere and HIF are en- 
gaged and this occurs at low densities. Thus the PC relation 
for SMC Cepheids should show a much reduced change in 
the PC slope in going from short to long period Cepheids 
and this is what we see in the observations: at maximum 
light, Galactic and LMC Cepheids display a PC(max) rela- 
tion with a slope that decreases to zero in going from short 
to long period (Paper I & III) but SMC Cepheids only show 
a reduced PC slope at maximum light. It may be possible 
to test this assertion spectroscopically and this will be the 
subject of a future investigation. Because the PC relation 
at mean light is an average of the PC relation at different 
phases, if there is a discontinuity or not at a certain phase, 
this effect will be transmitted to the mean light PC relation 
and hence to the mean light PL relation. 



4 CONCLUSION AND DISCUSSION 

In this paper we construct the full amplitude models ap- 
propriate for the SMC Cepheids to study the PC/ AC re- 
lations and the HIF-photosphere interaction. In short, the 
SMC models constructed in this paper do a reasonable job 
of matching the observed PC relations at minimum and 
mean light. The greatest discrepancies occur for short pe- 
riod Cepheids at maximum light. These SMC models tend 
to not be driven to higher temperatures and hence bluer 
(V — I) colours, though they do fall inside the outer enve- 
lope of the observed data points. 



From the previous papers in the series we note that the 
empirical LMC PC relation at minimum light has a statis- 
tically significant slope change or "break" at log(P) ~ 1.0, 
but has the smallest scatter over all phases. The LMC PC 
relation at maximum light has the largest scatter over all 
phases. It also shows a statistically significant change of 
slope across log(P) ~ 1.0 but has a flat slope (i.e. close 
to zero) when only the long period Cepheids are considered. 
The empirical Galactic PC relation displays a break at 10 
days with a flat relation thereafter at the maximum light 
and a linear relation at minimum light. The empirical SMC 
PC relation also display a marginal break at maximum light 
but is linear at minimum light. The PC relation is flatter for 
the SMC Cepheids with log(P) > 1 than the short period 
SMC Cepheids at the maximum light. 

In terms of the empirical AC relation, the AC rela- 
tions at maximum/minimum light are significantly broken 
(at log[P] = 1.0) for all three galaxies (see Paper I & III). 
At maximum light, higher amplitude Cepheids in all three 
galaxies are driven to hotter temperature and hence bluer 
colours for all period ranges (with negative slope in the AC 
relation). Nevertheless the slope of these empirical AC rela- 
tions at maximum light becomes shallower for the long pe- 
riod Cepheids. However, the opposite behavior is found at 
the minimum light as the long period Cepheids are driven 
to cooler temperature as the amplitude increases in all three 
galaxies. For short period Cepheids, the same behavior is 
found for the Galactic Cepheids but the LMC and SMC 
short period AC relations are nearly flat at minimum light. 
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Figu re 5. The pl ots of A as function of log(P). The open and squares and open circles are the models calculated with lBono et all (2000) 
and [chiosi (198sj) ML relation, respectively. The dashed lines represent (roughly) the outer boundary of the HIF. 
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Figure 6. Log of the photospheric density (defined as 1/V where V is the specific volume) plotted against period for the SMC models. 
The symbols are the same as in Figure [5] 



In all cases the behavior at mean light is intermediate be- 
tween the properties at maximum/minimum light. 

The Saha ionization equation, which is used in most 
models of Cepheid envelopes, is such that when the photo- 
sphere and HIF interact at low densities, the photospheric 
temperature will be almost independent of pulsation period. 
This is because when the HIF and photosphere are interact- 
ing or engaged, the photospheric temperature is essentially 
the same as the temperature at which hydrogen ionizes, and 
not dependent on the pulsation period. Hence the colour, 
either in (V — I) or in (B — V), is also almost independent 
of the pulsation period. Moreover, as the period increases, 
the L/M ratio increases which, in turn, changes the location 
of the HIF, in terms of the mass distribution, in the Cepheid 



envelope. This can change the phase at which the HIF and 
photosphere interact. 

One possible difference between Galactic, LMC and 
SMC models is the photospheric density at minimum light. 
For Galactic models, the L/M ratio causes the low density 
photosphere-HIF interaction to only occur at/around the 
maximum light, leading to a stronger relation between pe- 
riod and colour at minimum light. At other phases, the HIF 
and photosphere are well separated. For LMC models, the 
L/M ratio forces the HIF-photosphere interaction to occur 
at all phases for short period Cepheids at higher densities 
and at maximum light for longer period Cepheids at low 
densities. At phases close to the minimum light, the HIF- 
photosphere become disengaged for longer period Cepheids. 
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Table 6. Summary of various properties at the maximum and minimum light for the Galactic, LMC and SMC Cepheids and models. 
Amp and T refer to the amplitude and temperature, respectively. 



Galaxy Photospheric Density HIF-Photosphere Interaction Observed PC Relation 



Observed AC Relation 



GAL Low for log(P) > 0.7 



LMC 



SMC 



GAL 



LMC 



SMC 



Low for log(P) > 1.0 
High for log(P) < 1.0 

High for log(P) > 1.0 
High for log(P) < 1.0 



Yes for all periods 
Yes for all periods 
Yes for all periods 



Low for log(P) > 0.7 No for log(P) > 0.7 



Low for log(P) > 1.0 
High for log(P) < 1.0 

Low for log(P) > 1.5 
High for log(P) < 1.5 



No for log(P) > 1.0 
Yes for log(P) < 1.0 

Yes for all periods 



At Maximum Light 

Flat for log(P) > 0.8 with 
significant break at log(P) = 1.0 

Flat for log(P) > 1.0 with 
significant break at log(P) = 1.0 

^Flatter for log(P) > 1.0 with 
marginal break at log(P) = 1.0 

At Minimum Light 

Non-flat slope for all periods 
& no significant break 

Non-flat slope for all periods, 
significant break at log(P) = 1.0 

Non-flat slope for all periods 
& no significant break 



Higher Amp. — * hotter T. & 
gets shallower for log(P) > 1.0 

Higher Amp. — » hotter T. & 
gets shallower for log(P) > 1.0 

Higher Amp. hotter T. & 
gets shallower for log(P) > 1.0 



Higher Amp. — » cooler T. & 
gets shallower for log(P) > 1.0 

Flat for log(P) < 1.0 & 
gets steeper for log(P) > 1.0 

~Flat for all periods 



This leads to change in the slope of PC relation at phases 
around minimum light for the long and short period LMC 
Cepheids. For SMC models, the HIF-photosphere interac- 
tion occurs throughout the pulsation cycle, but only occurs 
at low density at phases close to, and including, the maxi- 
mum light for the long period models. This leads longer pe- 
riod SMC Cepheids to have a flatter PC relation slope than 
the short period SMC Cepheids at the maximum light. The 
high density photosphere-HIF interaction occurs at other 
phases and ensures that there is no significant change of 
slope in the PC relation as a function of period. The low 
density HIF-photosphere interaction at maximum light is 
the reason for the observed flat PC relation at maximum 
light for the Galactic and LMC Cepheids and the flatter 
PC (max) relation for longer period SMC Cepheids. In Ta- 
bleEl we summarize some of the properties discovered in this 
series of papers for the Galactic, LMC and SMC Cepheids 
and models. 

Because these changes depend on the L/M ratio and 
because the ML relation is dependent on metallicity, the 
above argument outlines a mechanism whereby the observed 
changes in the PC relation between Galactic, LMC and SMC 
Cepheids could occur. Moreover these changes do not have 
to be monotonic with metallicity. Since the mean light PC 
relation is the average of the PC relation at different phases, 
changes in the slope of the PC relation, can, in principle, 
affect the mean light PC relation (see Paper I to IV). 

In addition to our interest in understanding the physics 
behind the multi-phase PC and PL relations as a function 
of metallicity (in Galaxy, LMC and SMC), we are also inter- 
ested in the changes of light curves shape in different galax- 
ies. This is because the AC relations are essentially equiva- 
lent to plotting the light curve shape against temperature (or 
colour). Previous work has established that Principal Com- 
ponent Analysis (PCA) is a very efficient way of describing 
light curve structure with upwards of 80% of structure vari- 
ation being explained by just the first P rincipal Component 
l|Kanbur et al.ll2002l ; lTanvir et alj|2005l ). Since the first PCA 



coefficient correlates well with amplitude, therefore there 
will be a correlation between the am plitude and light curve 
structure. We remark in passing that iKeller fc Wood! (|2006T ) 
mention a correlation between amplitude and effective tem- 
perature in their hydrodynamic models of Cepheids. It can 
be seen from this work and previous papers in this series 
that a correlation between amplitude and effective temper- 
ature is due to a more extreme correlation at maximum or 
minimum light: the physics of which has been explained in 
this paper and earlier papers in the series. 
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